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Abstract 

In order to account for competition and interplay of localized 
and itinerant magnetic behaviour in correlated many body systems 
with complex spectra the various types of spin-fermion models have 
been considered in the context of the Irreducible Green's Functions 
(IGF) approach. Examples are generalized d — f model and Kondo- 
Heisenberg model. The calculations of the quasiparticle excitation 
spectra with damping for these models has been performed in the 
framework of the equation- of-motion method for two-time tempera- 
ture Green's Functions within a non-perturbative approach. A unified 
scheme for the construction of Generalized Mean Fields (elastic scat- 
tering corrections) and self-energy (inelastic scattering) in terms of 
the Dyson equation has been generalized in order to include the pres- 
ence of the two interacting subsystems of localized spins and itinerant 
electrons. A general procedure is given to obtain the quasiparticle 
damping in a self-consistent way. This approach gives the complete 
and compact description of quasiparticles and show the flexibility and 
richness of the generalized spin-fermion model concept. 



'International Journal of Modern Physics B13, N 20, (1999) 2573-2605 



^E-mail:kuzemsky@thsunl.jinr.ru ; http://thsunl.jinr. ru^ kuzemsky 



1 



1 Introduction 



The existence and properties of localized and itinerant magnetism in metals, 
oxides and alloys and their interplay is an interesting and not yet fully un- 
derstood problem of quantum theory of magnetism. The behaviour and the 
true nature of the electronic and spin states and their quasiparticle dynamics 
are of central importance to the understanding of the physics of correlated 
systems such as magnetism and Mott-Hubbard metal-insulator transition 
in metal and oxides, magnetism and heavy fermions (HF) in rare-earths 
compounds, high-temperature superconductivity (HTSC) in cuprates and 
anomalous transport properties in perovskite manganates. This class of sys- 
tems are characterized by the complex, many-branch spectra of elementary 
excitations and, moreover, the correlations effects are essential. 
Recently there has been considerable interest in identifying the microscopic 
origin of quasiparticle states in such systems and a number of model ap- 
proaches have been proposed. A principle importance of these studies is con- 
cerned with a fundamental problem of electronic solid state theory, namely 
with the tendency of 3d electrons in transition metal compounds and 4f 
electrons in rare-earth metal compounds to exhibit both localized and de- 
localized behaviour. The interesting electronic and magnetic properties of 
these substances are intimately related to this dual behaviour of electrons. 
In spite of experimental and theoretical achievements still it remains much 
to be understood concerning such systems. 

Many magnetic and electronic properties of these materials may be inter- 
preted in terms of combined spin-fermion models (SFM), which include the 
interacting spin and charge subsystems. This add the richness to physical 
behaviour and brings in significant and interesting physics, e.g. the bound 
states and magnetic polarons, HF, and colossal negative magnetoresistance. 
The problem of the adequate description of quasiparticle many-body dynam- 
ics of generalized spin-fermion models has been studied intensively during the 
last decades, especially in the context of magnetic and transport properies 
of rare-earth and transition metals and their compounds [1] - ||, magnetic 
semiconductors ||, ||, interplay of magnetism and HF J7|,[@] , HTSC |J 
- |14| and magnetic and transport properties of perovskite manganates ||T5fl , 



16] , [|J. Variety of metal- insulator transitions and correlated metals phe- 
nomena in <i(/)-electron systems as well as the relevant models have been 
comprehensively discussed recently in Ref. ||. 

The basic theory of the physical behaviour of SFMs has been studied mainly 
within mean field approximation. However many experimental investiga- 
tions call for a better understanding of the nature of solutions (especially 
magnetic) to the spin-fermion and related correlated models, such as t — J, 
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Kondo-Heisenberg model, etc fl"2|| . 



In the previous papers we set up the formalism of the method of the Irre- 
ducible Green's Functions (IGF) JT7| -||25|. This IGF method allows one to 
describe the quasiparticle spectra with damping for the systems with com- 
plex spectra and strong correlation in a very general and natural way. This 
scheme differs from the traditional method of decoupling of the infinite chain 
of the equations [26f and permits to construct the relevant dynamical so- 
lutions in a self-consistent way on the level of the Dyson equation without 
decoupling the chain of the equation of motion for the GFs. 
In this paper we apply the formalism to consider the quasiparticle spectra 
for the complex systems, consisting of a few interacting subsystems. It is the 
purpose of this paper to explore more fully the notion of Generalized Mean 
Fields (GMF) Jl2| which may arise in the system of interacting localized 
spins and lattice fermions to justify and understand the "nature" of relevant 
mean-fields and damping effects. 

It is worthy to emphasize that in order to understand quantitatively the 
electrical, thermal and superconducting properties of metals and their alloys 
one needs a proper description an electron-lattice interaction too |27]]- |29| . 



A systematic, self-consistent simultaneous treatment of the electron-electron 
and electron-phonon interaction plays an important role in recent studies of 
strongly correlated systems [23]|. The natural approach for the description 
of electron-lattice effects in such type of compounds is the Modified Tight- 
Binding Approximation (MTBA) pB| , |29| |. We shall consider here the effects 
of electron-lattice interaction within the spin-fermion model approach. 
The purpose of this paper is to extend the general analysis to obtain the 
quasiparticle spectra and their damping of the concrete model systems con- 
sisting of two or more interacting subsystems within various types of spin- 
fermion models to extend their applicability and show the effectiveness of 
IGF method. 



2 Irreducible Green's Functions Method 

In this paper we will use the IGF approach which allows one to describe 
completely the quasi-particle spectra with damping in a very natural way. 
The essence of our consideration of the dynamical properties of many-body 
system with complex spectra and strong interaction is related closely with 
the field theoretical approach and use the advantage of the Green's functions 
language and the Dyson equation. It is possible to say that our method tend 
to emphasize the fundamental and central role of the Dyson equation for the 
single-particle dynamics of the many-body systems at finite temperatures. 
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In this Section, we will discuss briefly this novel nonperturbative approach 
for the description of the many-body dynamics of many branch and strongly 
correlated systems. The considerable progress in studying the spectra of el- 
ementary excitations and thermodynamic properties of many-body systems 
has been for most part due to the development of the temperature dependent 
Green's Functions methods. We have developed a helpful reformulation of 
the two-time GFs method |26] which is especially adjusted [|17] for the corre- 



lated fermion systems on a lattice and systems with complex spectra 0,0. 
The similar method has been proposed in Ref. |3(J for Bose systems( an- 



harmonic phonons and spin dynamics of pure Heisenberg ferromagnet). The 
very important concept of the whole method are the Generalized Mean 
Fields. These GMFs have a complicated structure for the strongly corre- 
lated case and complex spectra and do not reduce to the functional of the 
mean densities of the electrons or spins, when we calculate excitations spec- 
tra at finite temperatures. 

To clarify the foregoing, let us consider the retarded GF of the form |26| 



G r =« A(t), B(t') »= -i9{t - t') < [A(t)B(t% >,r] = ±l. (1) 

As an introduction of the concept of IGFs let us describe the main ideas of 
this approach in a symbolic form. To calculate the retarded GF G(t — t') let 
us write down the equation of motion for it: 

uG(u) =< [A, A + ] v > + « [A, H}_\A + . (2) 



The essence of the method is as follows [^OJ. It is based on the notion of the 
"IRREDUCIBLE" parts of GFs (or the irreducible parts of the operators, 
out of which the GF is constructed) in term of which it is possible, without 
recourse to a truncation of the hierarchy of equations for the GFs, to write 
down the exact Dyson equation and to obtain an exact analytical representa- 
tion for the self-energy operator. By definition we introduce the irreducible 
part (ir) of the GF 

ir « [A,H]-\A + »=« [A, H]_ - zA\A + » . (3) 

The unknown constant z is defined by the condition (or constraint) 

<[[A,Ht,A + ] v >=0 (4) 

From the condition (4) one can find: 

<[[A,HUA+] V > _M t 

<[A,A+} V > M [) 
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Here M and Mi are the zeroth and first order moments of the spectral 
density. Therefore, irreducible GF are defined so that they cannot be reduced 
to the lower-order ones by any kind of decoupling. It is worthy to note that 
the irreducible correlation functions are well known in statistical mechanics. 
In the diagrammatic approach the irreducible vertices are defined as the 
graphs that do not contain inner parts connected by the G°-line. With the 
aid of the definition (3) these concepts are translated into the language of 
retarded and advanced GFs. This procedure extract all relevant (for the 
problem under consideration) mean field contributions and puts them into 
the generalized mean-field GF, which here are defined as 



(6) 

To calculate the IGF ir « [A, H]-(t),A + (t') » in (2), we have to write 
the equation of motion after differentiation with respect to the second time 
variable t'. The condition (4) removes the inhomogeneous term from this 
equation and is a very crucial point of the whole approach. If one introduces 
an irreducible part for the right-hand side operator as discussed above for 
the "left" operator, the equation of motion (2) can be exactly rewritten in 
the following form 

G = G° + G°PG°. (7) 
The scattering operator P is given by 

P=(M )- 1 ir «[A,H]_\[A + ,H]_» ir (Mo)' 1 . (8) 

The structure of the equation (7) enables us to determine the self-energy 
operator M, in complete analogy with the diagram technique 

P = M + MG°P. (9) 

From the definition (9) it follows that the self-energy operator M is defined 
as a proper (in diagrammatic language "connected") part of the scattering 
operator M = {P) p . As a result, we obtain the exact Dyson equation for the 
thermodynamic two-time Green's Functions: 

G = G° + G°MG, (10) 

which has a well known formal solution of the form 

G =[(G )- 1 - M}- 1 ; M = G 1 -G- 1 (11) 

Thus, by introducing irreducible parts of GF (or the irreducible parts of the 
operators, out of which the GF is constructed) the equation of motion (2) 
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for the GF can be exactly (but using constraint (4)) transformed into Dyson 
equation for the two-time thermal GF. This is very remarkable result, which 
deserves underlining, because of the traditional form of the GF method did 
not include this point. The projection operator technique has essentially 
the same philosophy, but with using the constraint (4) in our approach we 
emphasize the fundamental and central role of the Dyson equation for the 
calculation of the single-particle properties of the many-body systems. It is 
important to note, that for the retarded and advanced GFs the notion of 



the proper part is symbolic in nature |2(J • However, because of the identical 
form of the equations for the GFs for all three types (advanced, retarded 
and causal), we can convert in each stage of calculations to causal GFs and, 
thereby, confirm the substantiated nature of definition (9)! We therefore 
should speak of an analogue of the Dyson equation. Hereafter we will drop 
this stipulation, since it will not cause any misunderstanding. It should be 
emphasized that the scheme presented above give just an general idea of the 
IGF method. The specific method of introducing IGFs depends on the form 
of operator A, the type of the Hamiltonian and the conditions of the prob- 
lem. The general philosophy of the IGF method lies in the separation and 
identification of elastic scattering effects and inelastic ones. This last point 
is quite often underestimated and both effects are mixed. However, as far as 
the right definition of quasiparticle damping is concerned, the separation of 
elastic and inelastic scattering processes is believed to be crucially important 
for the many-body systems with complicated spectra and strong interaction. 
From a technical point of view the elastic (GMF) renormalizations can ex- 
hibit a quite non-trivial structure. To obtain this structure correctly, one 
must construct the full GF from the complete algebra of the relevant op- 
erators and develop a special projection procedure for higher-order GF in 
accordance with a given algebra. 

It is necessary to emphasize that there is an intimate connection between the 
adequate introduction of mean fields and internal symmetries of the Hamilto- 
nian. In many-body interacting systems, the symmetry is important in clas- 
sifying of the different phases and in understanding of the phase transitions 
between them. The problem of finding of the ferromagnetic and antiferromag- 
netic "symmetry broken" solutions of the correlated lattice fermion models 
within IGF method has been investigated in Ref. P5(| . A unified scheme for 
the construction of Generalized Mean Fields (elastic scattering corrections) 
and self-energy (inelastic scattering) in terms of the Dyson equation has been 
generalized in order to include the presence of the " source fields" . The " sym- 
metry broken" dynamical solutions of the Hubbard model, which correspond 
to various types of itinerant antiferromagnetism has been discussed. This 
approach complements previous studies of microscopic theory of Heisenberg 
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antiferromagnet jnj and clarifies the nature of the concepts of Neel sublat- 
tices for localized and itinerant antiferromagnetism and "spin-aligning fields" 
of correlated lattice fermions. 



3 Quasiparticle Dynamics of the d — f Model 

3.1 Generalized d — f model 

The concept of the s(d) — f model play an important role in the quantum 
theory of magnetism , P?l . In this section we consider the generalized d — f 
model, which describe the localized 4/(5/)-spins interacting with <i-like tight- 
binding itinerant electrons and take into consideration the electron-electron 



and electron-phonon interaction in the framework of MTBA f28fl,p9 
The total Hamiltonian of the model is given by 

H = H d + H d ^f + H d ^ ph + H P h (12) 

The Hamiltonian of tight-binding electrons is given by 

ij o" ia 



This is the Hubbard model j31]]. The band energy of Bloch electrons e(k) is 
defined as follows 

Ui = N- 1 <k) Bxp[ik(Ri - Rj)}, 



where N is the number of the lattice sites. For the tight-binding electrons in 
cubic lattice we use the standard expression for the dispersion 

e(k) = 2 tfia) cos(ka a ) (14) 

a 

, where a a denotes the lattice vectors in a simple lattice with inversion centre. 
The term H d _f describes the interaction of the total 4f(5f )-spins with the spin 
density of the itinerant electrons 

H d -f = E J<?& = -JN- 1 ' 2 E T\SZl4 a a k+q - a + z a S z _ q a + ka a k+q(J } (15) 

i kq & 

where sign factor z a is given by 

Zff = (+or— ) for a = (1 or |) 
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and 



_\S_ q if a = + 

- q ~\s^ Str = - 



In general the indirect exchange integral J strongly depends on the wave 
vectors J(k; k + q) having its maximum value at k = q = 0. We omit this 
dependence for the sake of brevity of notations only 



For the electron-phonon interaction we use the following Hamiltonian |28 



H d -p h = VU $, k + q)Q ?v at +qa a ka (16) 

kq 



where 



V v (k,k + q) = ^2t(a a )e°(q)[sma a k - sina a (fc - q)) (17) 



a 



here qo is the Slater coefficient |28[ originated in the exponential decrease of 
the wave functions of <i-electrons, N is the number of unit cells in the crystal 
and M is the ion mass. The e u (q) are the polarization vectors of the phonon 
modes. 

For the ion subsystem we have 

H Ph = \ E( p i^V + " 2 (q")QiQ g v) (is) 

^ qu 

where P qu and Q qu are the normal coordinates and oj{qv) are the acoustical 
phonon frequencies. Thus, as in Hubbard model the d- and s(p)-bands 
are replaced by one effective band in our d — f model. However, the s- 
electrons give rise to screening effects and are taken into effects by choosing 
proper values of U and J and the acoustical phonon frequencies. 

3.2 Spin Dynamics of the d — f Model 

In this section, to make the discussion more concrete and to illustrate the 
nature of spin excitations in the d — f model we consider the double-time 
thermal GF of localized spins p6| , which is defined as 



G+-(k;t-t') =« SUt),Sz k (t') »= -i9(t-t') < lS+(t),Sl k (t')}_ >= 

r+oo 

1/2tt / c^exp(-^)G + -(A;;c«/)(19) 



The next step is to write down the equation of motion for the GF. Our 
attention will be focused on spin dynamics of the model. To describe self- 
consistently the spin dynamics of the d—f model one should take into account 
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the full algebra of relevant operators of the suitable " spin modes" , which are 
appropriate when the goal is to describe self-consistently the quasiparticle 
spectra of two interacting subsystems. This relevant algebra should be de- 
scribed by the 'spinor' (j^j ("relevant degrees of freedom"), according to IGF 
strategy of Section 2. 

Once this has been done, we must introduce the generalized matrix GF of 
the form 

(<<!'£:>> <<S=:>>M ( * ;w) <2o) 

Here 

°k = a k1 a k+qi; a k =Y1 a kl a k+q] 

q q 

To explore the advantages of the IGF in the most full form, we shall do the 
calculations in the matrix form. 

To demonstrate the utility of the IGF method we consider the following 
steps in a more detail form. Differentiating the GF << S£(t)\B(t') » with 
respect to the first time, t, we find 

* « S t lB »„= { 2N - 1/2 < S » >} + ,2!) 

■T7 X! << ^^k-q a p1 a p+ql ~~ Sk-q( a pl a p+q1 ~ a pl a p+ql) \B »ui 

pq 

where 

MS} 

Let us introduce by definition irreducible (ir) operators as 

( S k-q) = S k-q~ < S 0> 5 k,q (22) 

{ a pi a p+qi) lT = °pt a p+<jT — ^ a p^\ a p^ > ^<?,o 

Using the definition of the irreducible parts the equation of motion (21) can 
be exactly transformed to the following form 

(w - JN-\n^ - ni )) « S£\B » u +2JN' 1 < S* >« a+\SZ k » u = (23) 







X! 2(Sfr_ q ) tr a^a p+q i — S^_ q (a^a p+q ^ — a^a p+q i) ir \B » u 
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Xfpl - fp+kl) 



where 

n<j = Yl < a t a i° >= S fir = I^( ex P(/ 3e (? (J )) + !) 
? q q 

To write down the equation of motion for the Fourier transform of the GF 
<< a k (t),B(t') » we need the following auxiliary equation of motion 

(u + e(p) - e(p - k) - 2JN- 1/2 < S% > -UN- 1 fa - nj) « a+ T a P+fci |5 » w +(24) 
^ _1 (/ P T " /p+fci « °t\B » u +JiV- 1 / 2 (/ pT - / p+fci « S+\B » u = 

} ~ JN 1 ^ 2 ^ << ^-r( a pTS+»*T^P+fc,g — a ql a p+kl^p,q+r) " \B » ljJ — 
) qr 

JAT-V2 J- « (^ r ) ir (a^a p+kl 5^ q+r + a+ a 9+ri 5 p+fej(? )|5 »« + 

qr 

UN 1 ^ << (ap T a+ +rT a,j|a p+r+fe j — ap +rT a^_ ri a 9 ja p+fc j) ir |5 >>uj + 

qr 

« (V U (<1, k + p- q)a^a k+p _ ql - V u (q,p)a+ +ql a k+pl )Q q 

uq 

Let us use now the following notations: 

P g 

B p = JN x l 2 ^2[S^ r (a^a q+r ^5 p+ k,q ~ a ^i a p+ki^p,q+r) ir — 

qr 

{S z _ r ) ir (a^a p+fe j.5 Pig+r + a^a q+r i5p + k,q)] + 
UN' 1 '^2(a+- [ a q ~ +r] a q] a p+r+kl - a+ +r] a+_ rl a ql a p+ki ) ir + 

qr 

J2(V v (q, k + p- q)a^a k+p - qi - V v (q,p)a+ +ql a k+pl ); 

uq 

n 1 =w-JN- 1 fa-n i ); Vt 2 = 2JN~ 1 < S^>; 

^ P ,k = + e(p) - e(p + fc) - A); 
A = 2J7V- 1 / 2 < 5 2 > -EWT 1 ^ - nj; 

x?(fc^) = iv- 1 E (/yffci " /pT) ; 

^ P ,fe 



In the matrix notations the full equation of motion can be summarized now 
in the following form 



(26) 
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where 

6 / n x n 2 \ f /j- 1 tv 1 /2q 2 o , 

Q= l-JAr-v^/ (i-^)J- M o -jvxfj' (27) 

* (p) = ( o .-I 

To calculate the higher-order GFs in (26), we will differentiate the r.h.s. of 
it with respect to the second-time variable (f). Combining both (the first- 
and second-time differentiated) equations of motion we get the " exact" ( no 
approximation have been made till now) "scattering" equation 

QG(k; cu) = / + £ Hp)P(p, qMqm+y 1 (28) 

pq 

This equation can be identically transformed to the standard form (7) 

G = G + G PG (29) 

Here we have introduced the generalized mean-field (GMF) GF Go according 
to the following definition 

Go = IT 1 / (30) 
The scattering operator P has the form 

P = i- 1 J2Hp)Hp,q)Hq)i- 1 (3i) 

pq 

Here we have used the obvious notation 

A „ , f«A ir \A ir » «A ir \B i a r »\ . . 

P(M;W)= [«^» «B-\E^») (32) 

The operators A and B q follow from A and B q by interchange f— >J,, fc — > —A; 
and 5+ -> -5-. 

As shown in Section 2, equation (29) can be transformed exactly into a Dyson 
equation (10) by means of the definition (9). Hence, the determination of 
the full GF G has been reduced to the determination of Go and M. 

3.3 Generalized Mean-Field GF 

From the definition (30) the GF matrix in generalized mean-field approxi- 
mation reads 

Go = R { shNrf -n lNx f) ^ 
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where 

R=(l- U X t f )^i + WN^xf 

The spectrum of quasiparticle excitations without damping follows from the 

poles of the generalized mean- field GF (33). 

Let us write down explicitly the first matrix element G^ 1 

«S + \S- »°= UN- 1 ' 2 < gg > 

uj-JN-\n 1 -n l ) + 2J 2 N-^<St>(l-U X o f )- l Xo f 

(34) 

This result can be considered as reasonable approximation for description of 
the dynamics of localized spins in heavy rare-earth metals like Gd. (c.f. M, 



The magnetic excitation spectrum following from the GF (34) consists of 
three branches - the acoustical spin wave, the optical spin wave and the 



Stoner continuum [27]. In the hydrodynamic limit, k — > 0, to — > the GF 



(34) can be written as 

where the acoustical spin wave energies are given by 

_ n , a _ 1/2 E,^ + ^K%)Mg + (2A)- 1 E g (/ gT - /g|)(%(g)) 2 
h ^- Uk ~ 2iVV2 < SS > +(n T - n 4 ) 

and 



(36) 



< 3 >=< ss > [i + a^^ ]- (37) 

For s.c. lattice the spin wave dispersion relation (36) becomes 

E{k) = < S* > +(n T - n;))" 1 (38) 

2t 2 a 2 

(— r— 53(/<?T ~ /alX** sin (^a) + fey sm(q y a) + k z sin(g 2 a)) 2 - 
a 

^XX-ZVr + /g|)(fex cosg x a + fe 2 , cosg^a + A; 2 cosg 2 a)) 
<? 

In GMF approximation the density of itinerant electrons ( and the band 
splitting A) can be evaluated by solving the equation 

na = E < 4^ka >= Y.[ Q M^<k)+UN- 1 n_ a -JN- 1 / 2 < S* > -e^+l] -1 
k k 

(39) 
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Hence, the stiffness constant D can be expressed by the parameters of the 
Hamiltonian (12). 

The spectrum of the Stoner excitations is given by p| 

oj k = e(k + q)- e(q) + A (40) 

If we consider the optical spin wave branch then by direct calculation one 
can easily show that 

E opt (k) = E° opt + D(UE opt /JA - l)k 2 

E° opt = J(ni-n l )+2J<S*> (41) 

From the equation (33) one also finds the GF of itinerant spin density in the 
generalized mean field approximation 

""■■''-'^^ips (42) 

3.4 Dyson Equation for d-f model 

The Dyson equation (10) for the generalized d — f model has the following 
form 

G(k; lu) = G (k; ^) + £ G (p; uj)M(pq; uj)G(q; u) (43) 

pq 

The mass operator 

M(pq;uj) = P {p \pq]uj) 

describes the inelastic (retarded) part of the electron-phonon, electron-spin 
and electron-electron interactions. To obtain workable expressions for matrix 
elements of the mass operator one should use the spectral theorem, inverse 
Fourier transformation and make relevant approximation in the time correla- 
tion functions. The elements of the mass operator matrix M are proportional 
to the higher-order GF of the following (conditional) form 

( (ir) « (5 + )a fc+pCT1 a+ +g(72 a 9(T2 |(5-)a+ +S(T3 a+ 4 a r+SCT4 

For the explicit approximate calculation of the elements of the mass operator 
it is convenient to write down the GFs in (44) in terms of correlation functions 
by using the well-known spectral theorem [^6| : 

( y S + )dk+p CTi ap +qcr2 a qi j 2 \(S )ai+scT Z a ra A a r+so i »( tr ' ) ' p ) = 

1 r+ oc dcu' f +co 
— / -(expC&y) + 1) / exp(— iui t)dt 

< (5-(0)4+^3(*)<W«,+^(0l(^ + )«^ 1 < +? . 2 « ?CT2 > {ir) ' P ) (44) 
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Let us first consider the GF << A\A » appearing in Mn. Further in- 
sight is gained if we select the suitable relevant "trial" approximation for the 
correlation function on the r.h.s. of (44). In this paper we show that the ear- 
lier formulations, based on the decoupling or/and on diagrammatic methods 
can be arrived at from our technique but in a self- consistent way. Clearly 
the choice of the relevant trial approximation for the correlation function in 
(44) can be done in a few ways. The suitable or relevant approximations 
follow from the concrete physical conditions of the problem under consider- 
ation. We consider here for illustration the contributions from charge and 
spin degrees of freedom by neglecting higher order contributions between the 
magnetic excitations and charge density fluctuations as we did in the theory 
of ferromagnetic [S] and antiferromagnetic J32|| , [Q semiconductors. For 
example, a reasonable and workable one may be the following approximation 
of two interacting modes M 

«^»**»j^s £ L-f*-„ (45) 

^~ kkxk^k-ik^a ^ 1 ^ 

Im « S k _ k4 \S_ k _ k2 » LUl Im << a k3fT ak 3 +k4,(r\ a k 1 a a ki+k2<T »ui 2 

(exp(/3(cj 1 + cu 2 )) + l 



F(ui,u 2 , ) 



(exp(M) - l)(exp(/3u; 2 ) - 1) 



On the diagrammatic language this approximate expression results from the 
neglecting of the vertex corrections. 

The system of equations (43) and (45) form a closed self-consistent system 
of equations. In principle, one may use on the r.h.s. of (45) any workable 
first iteration-step forms of the GFs and find a solution by repeated itera- 
tions. It is most convenient to choose as the first iteration step the following 
approximations: 

Im « S+_ k4 \SZ k _ k2 »^ 2irN-V 2 < S* > 6(U! - E(k + k 2 ))5 ki . k2 ; (46) 

Im << a k3a a ks+k4Cr \a kicr a kl+k2Cr >>oj 2 ~ 
n(fk 3 a - /fci CT )<K^2 + e(^3^") - e(^3 + k A a))5 k . iM+k2 5 klM+ki 

Then, using (46) in (45), one can get an explicit expression for the Mix 

^ 414 2J2 V- t 1 + N W k + " f*A W + N(E(k + q))f pa {l - 2/ 3 

«A\A» ~^2^ 



p+qaj 



N 2 f£ to- E(k + q)-e(pa)+e(p + qa) 

(47) 

where 

e(ka) = e(k) +U < n_ a >; N(E(k)) = [exp(j3E(k)) - l]" 1 (48) 
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The calculations of the matrix elements M± 2 , M 2 i and M 2 2 can be done in 
the same manner, but with additional initial approximation for phonon GF 

« QkAQt »~ - ^ 2 (H) _1 (49) 

It is transparent that the construction of the GF << B p \B q » will consist of 
contributions of the electron-phonon, electron-magnon and electron-electron 
inelastic scattering. 

« B p \B q »=« B p \B q »P h ' e + « B p \B q » m ' e + « B p \B q » e " e 
As a result we find the explicit expressions for the GFs in mass operator 

« B p \B q »^- e = W £ uj-\ru) 

z rv a=± 

/ [l + N(au(rv)) - / p+g+ri ]/ pT + N(au(rv))f p+q+rl (l - 2/ pT ) 
^ u — (au(rv) - e(p |) + e(p + A; + r J.)) 

((^(r,p + k))\ >p+k - V»(r,p)V»(r,p + k)5 q , p+k+r ) + 
[1 + N(au(ru)) - / p+fci ]/ p+rT + N(au(rv))f p+ki (l - 2/ p+rT ) 
- (auo(rv) -e(p + r|) + e(p + /c |)) 

((^(r,p)) 2 ^ p+fe - V(r,p)V 1 '(r,p + k)S q )) (50) 

The contribution from inelastic electron-magnon scattering is given by 

2 J 2 

« B p \B q »— = -— < S5 > £ 
[1 + 7V(£(r)) - / p+fc+rT ]/ pT + iV(£(r))/ p+fc+rT (l - 2/ pT ) 



( 



- (£(r) - e(p ]) + e(p + k + r ])) 

+ 

[1 + N(E(r)) - f p+kl ]f p+rl + N{E{r))f p+kl {\ - 2f p+k[ ) , 

w - ( E (r) - e(p + r j) + e(p + fc |)) ^ 9 ' P+fc 1 J 

The term due to the electron-electron inelastic scattering processes becomes 

U 2 



« B p \B q » e " e = A . 2 



(E 



(1 - / p+fci )(l - fr+sjUrJp+kX + /p+fcj/r+.lCl ~ /rl)(l ~ fp+s^ 

w - (e(p + k i) - e(p + s T) + e(r + s j) - e(r |)) 

(1 — /p+fc+sj)(l ~ fr-sl)frlfpf + /p+fc+s|/r-st(l ~ /r|)(l ~ /pt) l 

o; - (e(p + fc + n j) - t(p ]) + e(r - s |) - e(r |)) 



y^r (! ~ /gl)( 1 ~ fp+kQfrjfp+q-r^ + /g|/p+fcj(l ~ /rj)(l ~ fp+q-m^) , 

r " w-(e(p + A;J.) + e(gJ,)-e(rJ.)-e(p + g-rT)) 

(1 _ /?+r|)(l - fp-rl)fplfq-kl + fq+rlfp-r^ ~ fp\){^ ~ fq-k] 



u - (e(q + r |) + e(p - r t) - e(p T) - e(q - k |)) 
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'])<W (52) 



In the same way for off-diagonal contributions we find 

2J 2 



<< A\B q »- 



N 2 



( 



[1 + N(E(r)) - /,+„]/, _ tt + JV(£(r))/, 



< 55 > £ 

r 

1 - 2fq-k]) 



uj - (£(r) + e(g + r f) - e(g - fc T)) 



[1 + iV(E(r)) - / gi ]/ g+r - fci + iV(E(r))/ gi (l - 2/ g+r , fci ) 
- (E(r) - e(q + r - k j) + e(g |)) 

And we have << B p \A >>=<< Ajfip+fc >>. 



+ 

) 



(53) 



3.5 Self-Energy and Damping 

Finally we turn to the calculation of the damping. To find the damping 
of the quasiparticle states in the general case, one needs to find the matrix 
elements of the mass-operator in (43). Thus we have 



G11 
G21 



G12 
G22 



Gon 
C021 



G012 
G022 



Mn M12 
M21 M22 



(54) 



From this matrix equation we have 



J* 



M n = ^2«A\A»; 



M21 = S^f? M " 1 «b p \a»; 
Ml2 = ^? k?rl «^,»; 
M22 = ^rnv2 EK^?,^ 1 « R p\B q »; (55) 

With (54) and (55) the GF G becomes 

6= ,\ . pffi-"" -<*k (56) 

<iet(G -' - M) V -(^ " M21) fe^ft - M n ) / 
Let us estimate the damping of magnetic excitations. From (56) we find 

« St\SZ k » u = {G nyi_z {kiUj) ( 57 ) 
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Here the self-energy E is given by 



E(k,u) =M U + 



i - u x f 

(J7V-V2 _ M 12 )(JN-^ 2 - M 21 )Nxt f ((l - U X f) + M 22 7V X ^) _1 (58) 

Let us consider the damping of the acoustical magnons. Considering only 
the linear terms in the matrix elements of the mass operator in (58), we find 
for small k and u 

« S+\S- k » w « (59) 

k k uj- E(k) - 27V-V2 < S z > s ( fc ^ v ' 



where 

1 - f/xo 7 (1 " U X d o 
Then the spectral density of the spin-wave excitations will be given as 



E(fc, W ) « M n + (M 12 + M 21 ) - rr A ; 7 + ; A V M 22 ( 60) 



- -ImG ll (k,uu + ie) = --Im « S£\S_ k >>, 

7T 7T 



2iV~ 1 / 2 < S* > T{k, 



(w - E{k) - A(k, uj)) 2 + T 2 (k, uj) 



(61) 



where 



A(k,u) = 2N' 1/2 <S%> ReE(k,uj) 
T(k, uj) = 2N~ 1 ^ 2 < Sq > ImE(fc, uj + e) (62) 

describes the shift and the damping of the magnons, respectively. 

Finally we estimate the temperature dependence of F(k,u) due to the mass 

operator terms in (58). Considering the first contribution in (58) we get 



Mn = Im « A\A » w « J 2 < S* > £((1 + N(E(k + q)) - f p<T )f p+qt7 + 

pqa 

N(E(k + q))f pfT (l - 2f p+qa ))S(uj - E(k + p) + e(p + q) - e(p))(63) 
Using the standard relations 

iV(£(g))|^ = (exptfDq 2 ) - l)" 1 (64) 
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we find 

ImMn ~ J < S* > j^Zn J d 3 p J q q 2 dq J d(cos0) 



1 rPE„ 

dx ~ T (65) 



2(3 D Jo expx — 1 

The other contributions to T(k,uj) can be treated in the same way, where 
M12, M21 and electron- magnon contribution to M22 are proportional to T, 
too. For the electron-phonon contribution to M22 we find 

ImM^ = Im«B p \B, »f ~ — [ x 2 dx T 3 (66) 

(3 J exp x — 1 

Hence, the damping of the acoustical magnons at low temperatures can be 
written as 

r(fc, cu)| fe ,^ ~ Ti + r 2 T + r 3 T 3 (67) 

where the coefficients Tj, (i = 1, 2, 3) vanish for k = cu = 0, and furthermore 
for J = 0. 



3.6 Charge dynamics of d-f model 

To describe the quasiparticle charge dynamics or dynamics of carriers of the 
d — f model (12) we should consider the equation of motion for the GF of 
the form 

G ka =« a ka \a+ a » (68) 
Performing the first time differentiation of (68) we find 

(u - e{k))G ka = 1 + « a Uq-v a p-° a k+w\ a to » ~ 

pq 

JJT/2^2( << S-q a k+q-a\ata >> + Z a « S^ q a k+qa \a^ a >>) + 
q 

Y / V a (k-q,k) «a k „ qa Q qu \at (T » (69) 

qua 

Following the previous consideration we should introduce the irreducible GFs 
and perform the differentiation of the higher-order GFs on second time. Us- 
ing this approach the the equation of motion (69) can be exactly transformed 
into the Dyson equation 

G kc {u) = G°» + G ka (cu)M ka (cu)G kcT (cu) (70) 
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where 

e°(ka) = e(k) - z aJ ^ < S z > +^n_ CT (71) 

Here the mass operator has the following exact representation 

M k „{u) = M|» + Ml~ m {u) + M e k : ph (co) (72) 



where 



U 2 

M£{u)) = ^ J2( (ir) « 4+ q -^ a p- a p+^\ a t+s-^r-.ai_ sa » (4r) ' p )(73) 



N 2 

pqrs 



qs 

(W« 5i ? G fe+9(J |^a+ +S(T »^))(74) 
ACV) = E E *S(P - !.P)<(P.P + q)( (ir) « Q q „a p „ qu \Q s ,a + p+q(J »<*>*)(75) 

As previously, we express the GF in terms of the correlation functions. In 
order to calculate the mass operator self-consistently we shall use the "pair" 



approximation [17l,||22|| for the M ee and approximation of two interacting 
modes for M e ~ m and M e ~ ph 0, | fffi| . Then the corresponding expressions 
can be written as 

U 2 f duiduj 2 duj 3 



ly PQ J 



where 



and 



UJ + LO\ — L0 2 — ^3 

F ee (uj u uj 2 ,uj 3 ) 

gp +q -aMgk+p,aM9p,-aM (76) 



g ka (u) = —Im « a ka \a^ a » L 

7T 



F ee (u u u 2 ,u 3 ) = (/(<*) (i - f(uj 2 ) - fM) + fMfM) 

Let us consider now the spin-electron inelastic scattering. As previously, we 
shall neglect of the vertex corrections, i.e. correlation between the propaga- 
tions of the charge and spin excitations. Then we obtain from (74) 

J 2 f duj\duj 2 



= ^E / 12 F em (-^ 

JM q J uo — uoi — oo 2 



(g k+Pi - a {u 2 ){—Im « Sl q \S q a » Ul ) + g k+Pt<T {u 2 ){-^Im « S z q \S z _ q »^))(77) 
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where 

F em (Lu u u 2 ) = (l + N(uj 1 )-f(u; 2 )) 

And finally we shall find the similar expression for electron-phonon inelastic 
scattering contribution (75) 

Mt7 h {u) = £ - Ml 2 / F e ~ ph (ui, lo 2 ) 

^ J LO - LOi - UJ 2 

g p -q,a(ui)(—Im « QqvlQ+v » Ua ) (78) 
7T H 

where 

F e - ph (u 1 ,u J2 ) = (l + N(LU 2 )- fM) 

Equations (70), (76), (77) and (78) form a closed self-consistent system of 
equations for one-fermion GF of the carriers for a generalized spin-fermion 
model. To find explicit expressions for the mass operator (72) we choose for 
the first iteration step in (76) - (78) the following trial approximation 

g ka (u) = 6(u - e°(kcx)) (79) 

Then we find 

]lf ee (l A — ^ \ " /p+go-(l ~ fk+pg ~ fq-cr) "t~ fk+pcrfq-cr ( %Q) 

kA '~ N^LU + e (q,-a)-e°(p + q,a)-e (k + pa) 1 ' 

For the initial trial approximation for the spin GF we take the expression 
(46) in the following form 

—Im « S°\S ZZ >>~ z^N- 1 ' 2 < S* >)8(u - z a E{q)) (81) 

7T 

Then we obtain || 

M e- m(l A _ 2J 2 < gg > \- fk +g ,i + N(E(q)) 



^ w - e 0( k + q ,l)- E(q) ' 

lu-e-m/ \ 2J 2 <gg> l-f k ^ + N(E(q)) 
kl ()= \u-e\k-qA)-E{q) (§2) 

This result is written for the low temperature region, when one can drop the 
contributions from the dynamics of longitudinal (zz) GF which is essential 
at high temperatures and in some special cases. 

In order to calculate the electron-phonon term (78) we need to take as initial 
approximation the expressions (49) and (79). We then get 



M e~ P h ( \ _ \ - \Vujp- q,P)\ 2 f 1 - fk-g,a + N(u(qu)) /fc- g , g + N(u(qv)) \ 

ka {U) ^ 2u{qv) \uj-t\k-q,])-uo{qv) + u-t\k-q,])+u{qvy 
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where 



An°t 2 

K(p - q,p}\ 2 = £ -^(sma a p- sina a (p- q)fK{q)\ 2 (84) 



Then analysis of the electron-phonon term can be done as in Ref. p8 |. 

For the fully self-consistent solution of the problem the phonon GF can be 

easily calculated too. The final result is 

« QkAQt »= - uj 2 (kv) - n^M)- 1 (85) 

where the polarization operator II has the form 

IMu,) = £ \Vu{q~ k, q)\ 2 J q - k 'l~ \ U (86) 

The above expressions were derived in the self-consistent way for the gener- 
alized spin-fermion model and for finite temperatures. 

It is important to note that to investigate the spin and charge dynamics in 
doped manganite perovskites the scheme described above should be modified 
to take into account the strong Hund rule coupling in these systems but it 
deserve of separate consideration. In the present paper to show clearly the 
advantage of the IGF approach we shall consider another interesting example, 
the dynamics of carriers for the Kondo-Heisenberg model. 



4 Dynamics of Carriers in the Spin-Fermion 
Model. 

4.1 Hole Dynamics in Cuprates 

To show the specific behaviour of the carriers in the framework of spin- 
fermion model we shall consider a dynamics of holes in HTSC cuprates. A 
vast amount of theoretical searches for the relevant mechanism of high tem- 
perature superconductivity deals with the strongly correlated electron mod- 



els |L2| . Much attention has been devoted to the formulation of successful the- 
ory of the electrons (or holes) propagation in the Cu02 planes in copper ox- 
ides. In particular, much efforts have been done to describe self-consistently 



the hole propagation in the doped 2D quantum antiferromagnet |34| - j47| . 
The understanding of the true nature of the electronic states in HTSC are 
one of the central topics of the current experimental and theoretical efforts in 



the field [12], [40]. Theoretical description of strongly correlated fermions on 



two-dimensional lattices and the hole propagation in the antiferromagnetic 
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background still remains controversial. The role of quantum spin fluctua- 
tions was found to be quite crucial for the hole propagation. The essence of 
the problem is in the inherent interaction (and coexistence) between charge 
and spin degrees of freedom which are coupled in a self-consistent way. The 
propagating hole perturbs the antiferromagnetic background and move then 
together with the distorted underlying region. There were many attempts 
to describe adequately this motion. However, a definite proof of the fully 
adequate mechanism for the coherent propagation of the hole is still lacking. 
In this paper we will analyse the physics of the doped systems and the true 
nature of carriers in the 2D antiferromagnetic background from the many- 
body theory point of view. The dynamics of the charge degrees of freedom 
for the Cu02 planes in copper oxides will be described in the framework 
of the spin-fermion (Kondo-Heisenberg) model |J, [38], using the approach 
described in Section 3. 



4.2 Hubbard model and t-J model 

Before investigating the Kondo-Heisenberg model it is instructive to consider 
the t — J model very briefly. The t — J model is a special development of the 
spin-fermion model approach which reflect the specific of strongly correlated 
systems. To remind this let us consider first the Hubbard model [[Jl[ . 



The model Hamiltonian which is usually refered as to Hubbard Hamiltonian 
is given by 

H = t ij a ta a j^ + Kiani-a (87) 

ija ia 

For the strong coupling limit, when Coulomb integral U S> W, where W is the 
effective bandwidth, the Hubbard Hamiltonian is reduced in the low-energy 
sector to t-J model Hamiltonian of the form 

H = - n^K+a^l - + H.C.) + Jj^SiSj (88) 

ija ij 

This Hamiltonian play an important role in the theory of HTSC Let us 
consider the carrier motion. The hopping at half-filling is impossible and this 
model describe the planar Heisenberg antiferromagnet. The most interesting 
problem is the behaviour of this system when the doped holes are added. In 
the t — J model (U — > oo) doped holes can move only in the projected space, 
without producing doubly occupied configurations (< n-j- > + < nj. >< 1). 
There is then a strong competition between the kinetic energy of the doped 
carriers and the magnetic order present in the system. According to Ref. |37] . 
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it is possible to rewrite first term in (88) in the following form 



H t = t 



<ij> 



af^Si Sf a it + at St a,j + h.c 



hi 



39) 



This form show clearly the nature hole-spin correlated motion over antifer- 
romagnetic background. It is follows from (89) that to describe in a self- 
consistent way a correlated motion of a carrier one need to consider the 
following matrix GF Function: 



G(i,j) 



/ « a 

« Sr 

\«st 



<T » 



a Jl » 



j.T 



>> 



aJ T >> 
a+ » 



« at 

« a a\ a J[ » 
«S- 

« si 



\ a ii » 

Hi » 



« a^\Sf » 
« a a |5+ » 

« sr\sf » 

« st\st » 



« a,tf\Sj >> \ 
« a ii 



\S~» 
« Si \Sj » 



« Sl\Sj >> J 
(90) 

It may be shown after straightforward but tedious manipulations by using 
IGF method of Section 2 that the equation of motion (2) for the GF (90) 
can be rewritten Dyson equation (10) 

G(i,j;uj) = G (i,j;uj) + J2 G o(hm;uj)M(m,n;uj)G(n,j;uj) (91) 



The algebraic structure of the full GF in (91) which follows from (11) is rather 
complicated. For clarity, we illustrate some features by means of simplified 
problem. 



4.3 Hole Spectrum of t — J model 

It is well known [[5|],|4^] how to write down the special ansatz for fermionic 
operator as a composite operator of dressed hole operator and spin operator 
for the case J ^> t. The hole operator hi corresponding to fermion operator 
a+ on the spin- up sublattice using the ansatz = hiSf and similarly for 
spin-down sublattice have been introduced ( for a recent discussion see e.g. 
Ref. [0). Then the Hamiltonian (89) obtain the form 

H t = tY,I ii h+h i {bt + h) ( 92 ) 

ij 

Here 6j and are the boson operators, which results from the Holstein- 
Primakoff transformation of spins into bosons. Equation (92) is not conve- 
nient form because of its non-diagonal structure. Caution should be exercised 
because the new vacuum is a distorted Neel vacuum. 
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The equation of motion (2) and (3) for the hole GF can be written in the 
following form 

u « hj\h~l » —t^Ijn < B nj >« h n \hl »= 

n 

Sjk + tJ^IjnC «h n B nj \ht») (93) 

n 

Here B nj = (6+ + bj). The "mean-field" GF (6) is defined by 

^2(u5ij - tlij < Bji >)G (i, k; uj) = 5 jk (94) 

i 

Note, that "spin distortion" < B mn > does not depend on (R m — R n ). Then 
the Dyson equation (91) becomes 

G{g, k) = G (g, k) + J2 G (g, j)M(j } l)G(l, k) (95) 

n 

where self-energy operator is given by 

M(j,l) = t 2 Y J I J n{ ir « h n B nj \h+B lm » ir )Imi (96) 

mn 

The standard IGF-method's prescriptions for the approximate calculation of 
the self-energy ( c.f. Section 3.4) , can be written in the form 

M(j, I; u) = t 2 £ I jn I ml r du ld u 2 1 + NM - fM (97) 

mn J -co LO — 0J\ — Ll>2 

{-Im « B nj \B lm » Ul )(-ImG(n,m;u 2 )) 

71 71 

It is worthy to note that the "mass" operator (97) is proportional to t 2 
. The standard iterative self-consistent procedure of IGF approach for the 
calculation of mass operator encounter the need of choosing as a first iteration 
"trial" solution the non-diagonal initial spectral function ImGo] in another 
words, there are no reasonable "zero-order" approximation for dynamical 
behaviour. The initial hole GF in paper p6l was defined as 



G (j,k;uj) = ^— (98) 
u + xe 

which corresponds to static hole, without dispersion. In contrast, the ap- 
proximation for the magnon GF yield the momentum distribution of a free 
magnon gas. After integration in (97) , the mass operator is given by an 
expression quite similar to the one encountered in papers |3"5|, where the 
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Bogoliubov-de Gennes equations has been derived. It can be checked that 
the present set of equations (95) - (97) gives the finite temperature general- 
isation of the results |36|j . As we just mentioned, one of its main merits is 
that it enables one to see clearly the " composite" nature of the hole states in 
an antiferromagnetic background, but in the quasi-static limit. The recent 
analysis [[15[],||7j show that the difficulties of the consistent description of the 
coherent hole motion within t — J model are rather intrinsic properties of 
the model and of the very complicated many-body effects. From this point 
of view it will be instructive to reanalyse the less complicated model Hamil- 
tonian, in spite of the fact that its applicability has been determined as the 
less reliable. 



4.4 Kondo-Heisenberg Model 

As far as the CwCVplanes in the copper oxides are concerned, it was ar- 
gued P|,|38|| that a relatively reasonable workable model with which one can 
discuss the dynamical properties of charge and spin subsystems is the spin- 
fermion ( or Kondo-Heisenberg) model Q. This model allows for motion of 
doped holes and results from d-p model Hamiltonian. We consider the in- 
teracting hole-spin model for a copper-oxide planar system described by the 
Hamiltonian 

H = H t + H K + Hj (99) 
where H t is the doped hole Hamiltonian 

Ht = ~ E ( ta t a i° + H - C -) = E <k)ata ka (100) 

<ij>a ka 

where af a and CLicr £1X6 the creation and annihilation second quantized fermion 
operators, respectively for itinerant carriers with energy spectrum 

e(q) = -Ucos(l/2q x )cos(l/2q y ) = t lx {q). (101) 

The term Hj in (99) denotes Heisenberg superexchange Hamiltonian 

Hj = e Js m s n = ^ E J(v)s q s- q (102) 

<mn> q 

Here S n is the operator for a spin at copper site r n and J is the n.n. superex- 
change interaction 

J(q) = 2J[cos(q x ) + cos{q y )} = J l2 {q) (103) 
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Finally, the hole-spin (Kondo type) interaction Hk may be written as (for 
one doped hole) 

H K = J2KaA = N- 1 ^^K{q)[S^a£,a^^ + z.Si q a&at ¥v \ (104) 

i kq & 

This part of the Hamiltonian was written as the sum of a dynamic(or spin- 
flip) part and a static one. Here K is hole-spin interaction energy 

K(q) = K[cos(l/2q x ) + cos(l/2q y )} = K 73 (q) (105) 

We start in this paper with the one doped hole model (99), which is consid- 
ered to have captured the essential physics of the multi-band strongly corre- 
lated Hubbard model in the most interesting parameters regime t > J, \K\. 
We apply the IGF method to this 2D variant of the spin-fermion model. 
It will be shown that we are able to give a much more detailed and self- 
consistent description of the fermion and spin excitation spectra than in 
papers - [O], including the damping effects and finite lifetimes. 



For a recent discussion of the one-dimensional Kondo- Heisenberg model and 
the classification of the ground-state phases of this model in the context of 



a fixed-point strategy see Ref. [18 . 



4.5 Hole Dynamics in the Kondo-Heisenberg Model 

The two-time thermodynamic GFs to be studied here are given by 

G(ka, t - t') =« a ka (t), at{t ! ) »= -i6{t - t') < [a ka (t), at a (t')] + > 

(106) 

x +-(mn,t-t') =« S+(t),S-(t') »= -i6(t-t') < [5+(t), > 

(107) 

In order to evaluate the GFs (106) and (107) we need to use the suitable 
information about a ground state of the system. For the 2D spin 1/2 quan- 
tum antiferromagnet in a square lattice the calculation of the exact ground 
state is a very difficult problem. In this paper we assume the two-sublattice 
Neel ground state. To justify this choice one can suppose that there are 
well developed short-range order (c.f.Ref. j49|) or there are weak interlayer 
exchange interaction which stabilize this antiferromagnetic order. According 
to Neel model, the spin Hamiltonian (102) may be expressed as [19f|,|32l 



H J= E H jaP S m J nP (108) 

<mn> a ,(i 

Here (a,/?) = (a, b) are the sublattice indices. 

To calculate the electronic states induced by hole-doping in the spin- fermion 
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model approach we need to calculate the energies of a hole introduced in the 
Neel antiferromagnet. To be consistent with (108) and (90) we define the 
single-particle fermion GF as 

G(ka,u) = ( <K a f^ltl >> << a ^t a) Htl >:> ) (109) 
v ' \<< a b (ka)\a+(ka) » « a b [ka)\aj [ka) » J v ' 

Note, that the same fermion operators a a (ia), annihilates a fermion with 



spin a on the (a)-sublattice in the i-th unit cell has been used in paper |10 
The equation of motion for the elements of GF (109) are written as 

^2(Lu5 ay - e a/3 (/c)) « a,y(ka)\ap(ka) »= 5 af3 - « A(ka,a)\a^ » 

7 

(110) 

where 

A(ka, a) = AT-V2 £ K{p){SZ u pa a a {k + p - a) + z a S z _ pa a a (k + pa)) (111) 
p 

We make use of the IGF approach (see Section 2) to threat the equation of 
motion (110). It may be shown that equation (110) can be rewritten as the 
Dyson equation (10) 

G(ka, u) = G (ka, uj) + G (ka, u)M(ka, u)G(ka, u) (112) 

Here Go(ka,uj) = Q^ 1 describes the behaviour of the electronic subsystem in 
the Generalized Mean-Field(GMF) approximation . The Q matrix have the 
form 

where 

e a (ka) = e aa (k) - z a N~ l l 2 £ K(p) < S z pa > 5 pfi = e aa (k) - z a KS z (114) 

p 

S z = N- 1 ' 2 < S* a > 

is the renormalized band energy of the holes. 

The elements of the matrix GF Go{ka,uj) are found to be 

G {ka,u) = 7T^ + 77TT ( 115 ) 

uj — e + (ka) uj — e_(ka) 

, u(ka)v(ka) u(ka)v(ka) rbaf , , fUR , 

G {ka,uj) = — — — — = G [ka, uj) (116) 

uj — e + (ka) cu — e-{ka) 
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G {ka,uj) = — - H — — (117) 

uj — e + (ka) uj — e-(ka) 

where 

u\ka) = 1/2(1 - z a ^)-v\ka) = 1/2(1 + z a ^); (118) 

e±(fc<r) = ±R(k) = ((e ab (k) 2 + K 2 S 2 Z ) 1/2 (119) 

the simplest assumption is that each sublattice is s.c. and e aa {k) = 0(a = 
a,b). In spite that we have worked in the GFs formalism, our expressions 
(115) -(117) are in accordance with the results of the Bogoliubov (u,v)- 
transformation for fermions, but, of course, the present derivation is more 
general. 

The mass operator M in Dyson equation (112), which describes hole-magnon 
scattering processes, is given by as a "proper" part of the irreducible matrix 
GF of higher order 

yf (h(j ; _ , ^ « A{ka,a)\A+{k(T,a) »^ ^ « A{ka,a)\A + {ka,b) »(" : 



( ir ) « A(ka,b)\A+(ka : a) »^ ^ « A(ka,b)\A + (ka,b) »^ 

(120) 

To find the renormalization of the spectra e±(ka) and the damping of the 
quasiparticles it is necessary to determine the self-energy for each type of 
excitations. From the formal solution (11) one immediately obtain 

G±(ka) = {u- e±(ka) - E^A^uOr 1 (121) 

Here the self-energy operator is given by 

E ± (ka,u;) = A ± M aa ± A^M^ + M ba ) + A T M bb (122) 

where 

\v 2 {ka)J 
A\ = u(ka)v(ka) 

Equations (121) determines the quasiparticle spectrum with damping {uj = 
E — iV) for the hole in the AFM background. Contrary to the calculations 
of the hole GF in Section 4.3, the self-energy (122) is proportional to K 2 but 
not t 2 (c.f.eqn. (97)) 



1 + NM-fM 

UJ — LO\ — UJ 2 

(KfiQ, ^i)9ap{k + q-v,uj 2 ) + F%(q, Ui)g a p{k + q, uj 2 )) 



„ J-oo UJ — UJi — UJ<? 
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Here functions N(ou) and f{w) are Bose and Fermi distributions, respectively, 
and the following notations have been used for spectral intensities 

F%{q,u) = -Urn « S\ a \Sl qP (124) 
g a/3 (k<7,uj) = --Im « a a (ka)\at(ka) » LU ; i,j = (+, -, z). 
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The equations (123) and (112) forms the self-consistent set of equations for 
the determining of the GF (109). Coupled equations (123) and (112) can be 
solved analytically by suitable iteration procedure. In principle, we can use, 
in the right-hand side of (123) any workable first iteration step for of the 
relevant GFs and find a solution by repeated iteration. 



4.6 Dynamics of Spin Subsystem 



It will be useful to discuss briefly the dynamics of spin subsystem of the 
Kondo-Heisenberg model. When calculating the spin wave spectrum of this 



model we shall use the approach of Ref. (19| where the quasiparticle dynam- 
ics of the two-sublattice Heisenberg antiferromagnet has been studied within 
IGF method. The contribution of the conduction electrons to the energy 
and damping of the acoustic magnons in the antiferromagnetic semiconduc- 



tors within IGF scheme have been considered in Refs. [32], The main 
advantage of the approach of paper |L9| was the using of concept of "anoma- 
lous averages" (c.f. p5| ) fixing the relevant (Neel) vacuum and providing a 
possibility to determine properly the generalized mean fields. The functional 
structure of required GF has the following matrix form 



<< s kb \s_ ka » 



(125) 



Here the spin operators S ka ^ refer to the two sublattices (a, b). 
The equation of motion for GF (125) after introducing the irreducible parts 
has the form 



K 

iVV2 



7 

< S z a >« a+\B » w =< [S+ a , B]> + « Cl r a \B » w (126) 



where the following notations have been used 

'S 



B 



S 



—ka 
-kb. 



a = (a, b) 
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(Jj 



a = 2(< sz > J + iv~ 1/2 £ J,K) = -<"S ; 

f = 2(< S 6 2 > J fc + iV- 1 / 2 £ J*- ff O = -cuf ; 



q 

A ab _ 2 < (^-qa)"'^^)"' > /-, 97 x 

^ " 2JVV2<5*> 10 

The construction of the irreducible GF << Cj! a \B » is related with the 
operators 

/nrir Air . pir 

9 

/lir- V 7fQ+fQ 2 V r Q+ / qz V r V r ■ 

^fcg — A7-1/2 2- ^iWqbWk-qa) ^k-qaK^qb) ) 1 

q 

K K 

B % ka = ~~/\/T/2 ^2(Sk-q a y ra p1 a p+ql + ^T7 XI Z <jSk-qa( a t<j a P+<lo) lr ( 12 ^) 
pq pqa 

With the aid of (24) and (25) the equation of motion for the mixed GF can 
be written as 

KN 1 / 2 

« a+\B »= — — xgfau)^ « Str\B » + 

Z 7 

K E- L «W B » (!29) 

p U p,k 



2N 1 / 2 

Combining the equations of motion (126) and (129) we find 

n s x(kuj) = i + £>i (130) 

where 



2S 2 
-25, 



(131) 



Then equation (130) can be transformed exactly into the Dyson equation for 
the spin subsystem 

X(ku) = x\kuj) + X \ku)M s (ku)x(ku) (132) 



29 



Here 

x°(ku) = n; 1 i (133) 

The mass operator of the spin excitations is given by the expression 



M s (ku) 



«c£i(c£) fr » «cr a \(c+r» 

asi V « c£ b \(c+ a r » « cn\(ct b y r » 



(134) 



We are interesting here in the calculation of the spin excitation spectrum 
in the generalized mean field approximation. This spectrum is given by the 
poles of the GF x° 

dettt s (ku) = (135) 

Depending of the interrelation of the parameters this spectra have different 
forms. For the standard condition 2t ^> KS Z we obtain for the magnon 
energy [R3] 



101 



±u k = ± oWl - 72(^) 2 =F 



K 2 S Z 



Xo{k,u k ) t 



^ ^m) ( 136 ) 



The acoustic magnon dispersion law for the k — > is given by 



± 



±D(T)\k\ 



where the stiffness constant |33 

1 



D(T) = zJSjl 



NS 



K 2 S Z 



z q 



AN k^O 



(137) 



:i38i 



The detailed consideration of the spin quasiparticle damping will be done in 
separate publication. Here we now proceed with calculating the damping of 
the hole quasiparticles. 



4.7 Damping of Hole Quasiparticles 

It is most convenient to choose as the first iteration step in (123) the simplest 
two-pole expressions, corresponding to the GF structure for a mean field, in 
the following form 

g a/3 (ka,uj) = Z + 5(u - t+(ka)) + Z^5(uj - t-{ka)) (139) 

where Z± are the certain coefficients depending on u(ka) and v(ka). The 
magnetic excitation spectrum corresponds to the frequency poles of the GFs 
(107). Using the results of Section 4.5 on spin dynamics of the present model, 
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we shall content ourselves here with the simplest initial approximation for 
the spin GF occurring in (123) 



2z a S, 



■F°s a (q, uj) = L + 5(u - z a u q ) - L-5(u + z a uj q ) (140) 



Here u q is the energy of the antiferromagnetic magnons (136) and L± are the 
certain coefficients (see |B|). We are now in a position to find an explicit 
solution of coupled equations obtained so far. This is achieved by using (139) 
and (140) in the right- hand-side of (123). Then the hole self-energy in 2D 
quantum antiferromagnet for the low-energy quasiparticle band t-(ka) is 

= g& Ey?( i + -/(*-(*-«)) + NM + nuk + v)) 

V ' ; 2N V u-u q -t-(k-q) u + uj q -t^(k + q) V ; 

+ 2_K*Sl y2 iVK +P )(l + N(u q )) + f(Uk + p))(N(u q ) - N(u q+P )) 
N 2 u + u q+p - u q - t_(k + p) 

Here we have used the notations 

Y 2 = (U q + V q ) 2 ; Y 2 = (U q U q+p - V q V q+p ) 2 

where the coefficients U q and V q appears as a results of explicit calculation 



of the mean- field magnon GF |19[, |33|] 



A remarkable feature of this result is that our expression (141) accounts for 
the hole-magnon inelastic scattering processes with the participation of one 
or two magnons. 

The self-energy representation in a self-consistent form (123) provide a pos- 
sibility to model the relevant spin dynamics by selecting spin-diagonal or 
spin-off-diagonal coupling as a dominating or having different characteristic 
frequency scales. As a workable pattern, we consider now the static trial 
approximation for the correlation functions of the magnon subsystem fl9| in 
the expression (123). Then the following expression is readily obtained 

M&ka, uj) = —J2 ;(< Sl° p S° qa > g aP (k + q - a, u') (142) 

iv q J—od OJ — OJ 

+ <(S z _ q(s y r (S z q J r >g Q e(k + qa 1 uj')) 
Taking into account (141) we find the following approximate form 

E (k<r,u) » 7T1T7Z^ m Z S I 1 -73(g)) 143 

It should be noted, however, that in order to make this kind of study valuable 
as one of the directions to studying the mechanism of HTSC the binding of 
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quasiparticles must be taking into account. This very important problem 
deserves the separate consideration. In spite of formal analogy of the our 
model (99) with that of a Kondo lattice, the physics are different. There is a 
dense system of spins interacting with a smaller concentration of holes. This 
question is in close relation with the right definition of the magnon vacuum 
for the case when K ^ 0. 

In this Section we has considered the simplest possibility, assuming that dis- 
persion relation e aa (k) = (a = a, b). In paper [^TJ a model of hole carriers 



in an antiferromagnetic background has been discussed, which explains many 
specific properties of cuprates. The effect of strong correlations is contained 
in the dispersion relation of the holes. The main assumption is that the influ- 
ence of antiferromagnetism and strong correlations is contained in the special 
dispersion relation, e(k), which was obtained using a numerical method. The 
best fit corresponds to[[4"T| 

e(Jfe) =-1.255 + 0.34 cos ^ cos ife tf + 0.13(cos2A; x + cos 2jfe,,) (144) 

As a result, the main effective contribution to e(k) arises from hole hopping 
between sites belonging to the same sublattice, to avoid distorting the anti- 
ferromagnetic background. 

Our IGF method is essentially self-consistent, i.e. do not depends on the spe- 
cial initial form for the hole propagator. For the self-consistent calculation 
by iteration of the self-energy (123) we can take as the fist iteration step the 
expression (139) with the dispersion relation (144) or another suitable form. 
This must be done for the calculation mean-field GF (113) and dispersion 
relation (119) too. 

To summarize, in Section 4 we have presented calculations for normal phase 
of HTSC, which are describable in terms of the spin-fermion model. We have 
characterized the true quasiparticle nature of the carriers and the role of 
magnetic correlations. It was shown that the physics of spin-fermion model 
can be understood in terms of competition between antiferromagnetic order 
on the CwC^-plane preferred by superexchange J and the itinerant motion 
of carriers. In the light of this situation it is clearly of interest to explore 
in details how the hole motion influence the antiferromagnetic background. 
Considering that the carrier-doping results in the HTSC for the realistic 
parameters range t 3> J, K, corresponding the situation in oxide supercon- 
ductors, the careful examination of the collective behaviour of the carriers 
for a moderately doped system must be performed. It seems that this be- 
haviour can be very different from that of single hole case. The problem 
of the coexistence of the suitable Fermi-surface of mobile fermions and the 
antiferromagnetic long range or short range order has to be clarified. 
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5 Conclusions 



We have been concerned in this paper with establishing the essence of quasi- 
particle excitations of charge and spin degrees of freedom within a generalized 
spin-fermion model, rather than with their detailed properties. We have con- 
sidered the generalized d — f model and Kondo-Heisenberg model as the most 
typical examples but the similar calculation can be performed for other anal- 
ogous models. To summarize, we therefore reanalyzed within IGF approach 
the quasiparticle many body dynamics of the generalized spin-fermion model 
in a way which provides us with an effective and workable scheme for con- 
sideration of the quasiparticle spectra and their damping for the correlated 
systems with complex spectra. The calculated temperature behaviour of the 
damping of acoustical magnons (67) can be useful for analysis of the expere- 
mental results for heavy rare-earth metals like Gd P[. The present analysis 
of the 2D Kondo-Heisenberg model complements the previous analytical [13[ 
and numerical studies, showing clearly the important role of the damping 
effects. 

We have considered from a general point of view the family of solutions for 
itinerant lattice fermions and localized spins on a lattice, identifing the type of 
ordered states and then derived explicitly the functional of generalized mean 
fields and the self-consistent set of equations which describe the quasiparticle 
spectra and their damping in the most general way. While such generality 
is not so obvious in all applications, it is highly desirable in treatments of 
such complicated problems as the competition and interplay of antiferromag- 
netism and superconductivity, heavy fermions and antiferromagnetism etc., 
because of the non-trivial character of coupled equations which occur there. 
The problem of the coexistence of HF and magnetism is extremely nontriv- 
ial 0,0 many-body problem and have no appropriate solution in spite that 
there are many experimental evidences of the competition and interplay of 
HF and antiferromagnetism ||. Both these problems are subject of current 
but independent research. 

Another development of the present approach is the consideration of the 
competition and interplay of itinerant and localized magnetism and anti- 
ferromagnetism of the doped manganite perovskites where the interrelation 
between parameters of the spin-fermion model is quite different and the new 
scheme of approximation should be invented. Especially, the situation , when 
Hund rule interaction is very large but finite should be carefully analyzed. It 
would be interesting to understand on a deeper level the relationship between 
different possible phase states in manganates and various ordered magnetic 
states within the generalized spin-fermion model. 

In conclusion, we have demonstrated that the Irreducible Green's Functions 
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approach is a workable and efficient scheme for the consistent description of 
the quasiparticle dynamics of complicated many body models. 
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